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Abstract 

We study the cosmology of supersymmetric models in which the supersymmetry 
breaking effects are mediated by gauge interactions at about the 10^ GeV scale. We 
first point out that the gravitino is likely to overdose the Universe in this class of 
models. This requires an entropy production, which prefers a baryogenesis mechanism 
at a relatively low temperature. The Affleck-Dine mechanism for baryogenesis is one of 
the possibilities to generate enough baryon asymmetry, but the analysis is non-trivial 
since the shape of the potential for the flat direction differs substantially from the 
conventional hidden sector case. To see this, we first perform a 2-loop calculation to 
determine the shape of the potential. By combining the potential with the supergrav¬ 
ity contribution, we then find that the Affleck-Dine baryogenesis works efficiently to 
generate sufficient baryon asymmetry. On the other hand, we also point out that string 
moduli fields, if present, are stable and their coherent oscillations overdose the Uni¬ 
verse by more than 15 orders of magnitude. One needs a very late inflationary period 
with an e-folding of A ^ 5 and an energy density of ^(10^ GeV)^. A thermal inflation 
is enough for this purpose. Fortunately, the Affleck-Dine baryogenesis is so efficient 
that enough baryon asymmetry can survive the late inflation. 


*This work was supported in part by the Director, Office of Energy Research, Office of High Energy 
and Nuclear Physics, Division of High Energy Physics of the U.S. Department of Energy under Contract 
DE-AC03-76SF00098 and in part by the National Science Foundation under grant PHY-95-14797. HM was 
also supported by the Alfred P. Sloan Foundation. AdG was supported by CNPq (Brazil). 





Disclaimer 


This document was prepared as an account of work sponsored by the United States Government. While 
this document is believed to contain correct information, neither the United States Government nor any agency 
thereof, nor The Regents of the University of California, nor any of their employees, makes any warranty, express 
or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any 
information, apparatus, product, or process disclosed, or represents that its use would not infringe privately 
owned rights. Reference herein to any specific commercial products process, or service by its trade name, trade¬ 
mark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or 
favoring by the United States Government or any agency thereof, or The Regents of the University of California. 
The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States 
Government or any agency thereof, or The Regents of the University of California. 


Lawrence Berkeley Laboratory is an equal opportunity employer. 



1 Introduction 


The promise of low-energy supersymmetry (SUSY) is to stabilize the hierarchy between the 
weak scale and a higher scale of new physics, e.g. the Planck scale (see, e.g., ID). Super- 
symmetry, however, has to be spontaneously broken because we do not see degenerate pairs 
of particles and their superpartners. Moreover, there are stringent phenomenological con¬ 
straints on the spectrum of superparticles such as the degeneracy among squarks or sleptons 
at the percent level (see, e.g., p[). Therefore, constructing viable mechanisms of supersym¬ 
metry breaking has been regarded as one of the most important issues in supersymmetry 
model building. 

Recently, the idea of generating supersymmetry-breaking masses via gauge interactions 
has attracted interests (Low-Energy Gauge Mediation, or LEGM) ii.a. In this scheme, 
the supersymmetry breaking effects appear in the supersymmetric Standard Model in the 
following manner. There is a sector which breaks supersymmetry dynamically at around a 
10^ GeV scale, and it generates supersymmetry breaking effects in the so-called messenger 
sector at around 10® GeV, which further induce supersymmetry breaking masses of order 
10® GeV in the supersymmetric Standard Model via ordinary gauge interactions. This 
mechanism guarantees the required degeneracy among squarks, sleptons at a high degree, 
and also generates the masses for scalars and gauginos at comparable magnitudes as desired 
phenomenologically. There have been active studies on the phenomenology of such models 
1^, 1^. On the other hand, there has been little discussion on cosmological consequences of 
this mechanism, except issues concerning stable particles in the messenger sector pj. Since 
the scheme completely differs from the conventional hidden sector scenario at high energies, 
early cosmology is expected to differ substantially as well. 

There are (at least) two ingredients in the LEGM models which may lead to a cosmology 
different from the hidden sector case. The first is a very small gravitino mass. Since super- 
symmetry is broken at around 10^ GeV, compared to around 10^° GeV in the hidden sector 
scenario, the gravitino mass is much lighter: 7713/2 ~ 100 keV compared to 100 GeV. The 
second is that the supersymmetry breaking effects “shut off” at high energies.^ In particular, 
the flat directions in the supersymmetric Standard Model have very different potentials at 
large field amplitudes. 

In this paper, we study the implications of the LEGM models to cosmology. In Section ^ 
we hrst discuss the cosmological constraints on light gravitinos mainly based on the analysis 
by two of the authors (TM and HM) and Yamaguchi |^. Then, in Section we pay 


particular attention to the estimate of the gravitino mass in the LEGM models, and argue 
that it is highly unlikely to be lighter than 2h^keV as required by cosmology. This point 
implies that there must be a substantial entropy production, which casts a concern on the 
baryon asymmetry. Therefore, we turn our attention to a possible mechanism of baryogenesis 
at a relatively low temperature, using the idea by Affleck and Dine |^. The important point 
in the LEGM models is that the SUSY breaking effects due to the messenger interaction “shut 
off” at high energies. Therefore, in Section 0, we performed a 2-loop calculation to determine 


*A similar effect was discussed in H in the context of the sliding singlet mechanism. 
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the shape of the potential for the flat direction. In Section we estimate the possible 
value of the baryon-to-entropy ratio which can be induced by the Affleck-Dine baryogenesis. 
For a sufficiently large amplitude of the flat direction, the potential is dominated by the 
supergravity contribution rather than the LEGM contribution, and we will see that the 
Affleck-Dine baryogenesis works well enough to explain the present value of the baryon 
asymmetry. Furthermore, in Section ^ we point out that the string moduli, if present 
within the LEGM models, cause a serious problem because they are stable and their coherent 
oscillations grossly overdose the Universe. However we also point out a possible solution to 
the problem. Since the Affleck-Dine baryogenesis is so efficient, the baryon asymmetry can 
survive the enormous entropy production required to dilute the moduli fields, possibly by 
thermal inflation |^. Finally, in Section]^, we summarize our conclusions. 


2 Cosmology of a Light Gravitino 


In this section, we briefly review the cosmology with a light stable gravitino If a stable 

gravitino is thermalized in the early Universe, and if it is not diluted by some mechanism 
(such as a late inflation and/or a substantial entropy production), its mass density may 
exceed the closure limit: D 3/2 < 1- Since the number density of the gravitino is fixed once 
the gravitinos are thermalized, the above argument sets an upper bound on the gravitino 
mass 


7713/2 < 2 h^keV : without dilution, ( 2 . 1 ) 

where h is the Hubble constant in units of lOOkm/sec/Mpc. In other words, if the gravitino 
mass is heavier than 2h^ keV, we need some mechanism to dilute the gravitino in order not to 
overdose the Universe. Since the gravitinos are produced more at a higher temperature, we 
obtain an upper bound on the maximal temperature, Tmax, from which the ordinary radiation 
dominated Universe starts. For example, in the inflationary Universe, Tmax corresponds to 
the so-called reheating temperature Trh which is typically higher than Trh ~ 10® GeV, if 
there is no significant entropy production after reheating. If Tmax turns out to be less than 
10® GeV or so, we judge that one needs a substantial entropy production below Trh- It is 
worth to recall that the recent measurements prefer h ~ 0.7 and hence the upper bound is 
about 1 keV. 

The crucial point about the light gravitino is that the interaction of the (longitudinal 
component of) gravitino becomes stronger as the gravitino mass gets lighter. This is because 
the longitudinal component of the gravitino behaves like the goldstino, whose interaction is 
proportional to {F)~^ ~ (r? 73 / 2 M*)“^, where M* = 2.4 x 10^® GeV is the reduced Planck 
scale. For the light gravitino, the interaction of the longitudinal component of the gravitino 
(~ the goldstino) xjj to the chiral multiplet (0, x) and to the gauge multiplet (d^. A) is given 

*In this paper, we assume the absolute conservation of i?-parity. 
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■(V’Xl) 0* + h.c., 


( 2 . 2 ) 


where m^, m^, and m\ represent the masses of 0, y, and A.[] As indicated in Eq. (^), the 
interaction of ip becomes stronger as the gravitino mass gets smaller. 

In the thermal bath, two types of the processes may contribnte to overprodnce the 
gravitino: one is the decay of the sparticle X into its superpartner X and the gravitino, 
X ip + X, and the other is the scattering processes, x + y ^ ip + z, where x, y, z are 
relevant (s)particles. The decay process is significant especially for the case 7713/2 ^ lOOkeV. 
The partial decay width of a sparticle X into the gravitino is estimated as 


T{X ^iP + X) 




4877 ml^2^1* ’ 


(2.3) 


with mj^ being the mass of X, and it becomes large as the gravitino mass gets small. This 
decay process produces gravitinos as 773/2 + 8 ^^ 773/2 = r(X ip + X)nx where H is the 
expansion rate of the Universe at the given time.f] Here and below, 773/2 is the number 
density of gravitinos in the Universe at a given time. If the gravitino mass is in the range 
2h^keV ^ 7773/2 ^ lOOkeV, the decay rate becomes so large that the decay process overpro¬ 


duces the gravitino once the sparticles are thermalized 
in this range, Tmax should be lower than about m- 


0 




Thus, if the gravitino mass is 
100 GeV - 1 TeV depending on the 
mass spectrum of superparticles, or the Universe is overdosed.^ 

If the gravitino mass is heavier than (9(100keV), the decay process becomes unimportant 
and the most important production mechanisms of gravitinos are scattering processes. In 
this case, the Boltzmann equation for the number density of the gravitino, 773 / 2 , is given by 


^•3/2 + 31/773/2 — Stot77rad) 


(2.4) 


where H is the expansion rate of the Universe, Stot is the thermally averaged total cross 
section, and rirad = (C(3)/7r^)T^. At high energies the hrst term in Eq. (|2.2|) becomes more 
significant than the second one, and hence Stot is as large as 0{glmQ^/my2^*)■ After a 
detailed calculation, we obtain 0 


E 


tot 


rsj 


5.9 




(2.5) 


^Here, tp represents the spin ^ field, though the gravitino has spin |. In the high energy limit, tp is related 
to the longitudinal (helicity ±1/2) component of the gravitino, bf /21 s-s iPi /2 ~ lm^/ 2 - 

Hf the gravitino number density is large, there is also a damping term because of the detailed balance, 
(r.h.s.) = r(A %p + — 773 / 2 / 773 ^ 2)1 where is the thermal equilibrium value of 773 / 2 - 

§A similar argument can be applied to the decay process of the particles in the messenger sector or the 
SUSY breaking sector. In that case, the decay rate becomes much larger since the parent particle is much 
heavier. Thus, if the particles in those sectors are thermalized, the lower bound on the gravitino mass 
becomes more stringent than ^ lOOkeV. 
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where gs and mcs are the gauge coupling constant and the gaugino mass for SU(3)c. Solving 
Eq. (|2.4|) , and taking account of the effect of the dilution factor, g^iT^nx)/g*{T) (where g*(T) 
is the number of the relativistic degrees of freedom in the thermal bath with temperature 
T), the number density of the gravitino is given by 


nz/2{T) 

^rad 




5'*(^max) H T=Tn] 

3 X 10-2 3*{T) / 7713/2 fniGs y / 

g*(Tinax) VlOOkeV/ VlTeV/ V 


lOTeVy 


( 2 . 6 ) 


Using g^{T^ IMeV) ~ 3.9, and 5 '*(Tmax) ~ 200, we obtain 

_ “ 3 /.%/. ^ J ^ (r^y (T^\ 

' Pc UOOkeV/ VlTeV/ UOTeV/ 


(2.7) 


and the condition ^ 3/2 < 1 sets an upper bound on Tmax- In summary, the upper bound on 
Umax is given by 


T < 

-‘max ^ 


lOOGeV 


ITeV 


: 2/i^keV ^ 1713/2 ^ lOOkeV 


lOTeV X *3 > lookeV 


( 2 . 8 ) 


The above constraints are summarized in the Fig. |l[ As one can see, the upper bound on 
Umax is much lower than the usual reheating temperature after ordinary inflation, Trh ~ 10 ®. 
To reduce the number density of the gravitino, therefore, a large entropy production is 
required. 


3 Light Gravitino in the LEGM Models 

We discussed cosmological constraints on a light stable gravitino in the previous section, 
and showed that one needs to dilute gravitinos produced in the early Universe somehow if 
1713/2 ~ 2h.^ keV. In this section we estimate the gravitino mass in the LEGM models carefully 
and hnd it unlikely to be below 2h^ keV. 

In the scheme of the LEGM models, there are three sectors: the dynamical SUSY breaking 
(DSB) sector which originally breaks SUSY, the ordinary sector which consists of the particles 
in the minimal SUSY standard model (MSSM), and the messenger sector which mediates the 
SUSY breaking from the DSB sector into the ordinary sector. The scales for these sectors 
have a large hierarchy, since they are related by loop factors: Adsb ^ Amess S> Msusy, 
where Adsb, A^ess and Msusy represent the scales for the dynamical SUSY breaking sector, 
messenger sector, and the ordinary sector (~ electroweak scale), respectively. 

In the LEGM models a gauge interaction, which becomes strong at the scale Adsb, 
induces a non-perturbative superpotential. Due to non-perturbative effects, F-components 
of chiral multiplets in the SUSY breaking sector acquire non-vanishing VEVs, (Fq) ^DSB, 
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Figure 1: The upper bound on T^ax as a function of the gravitino mass from the requirement 
that the relic stable gravitinos do not overdose the Universe. We take the Hubble param¬ 
eter to be Hq = 100 Mpc/km/sec. There is no constraint below m ^/2 = 2 keV, which is 
represented by the vertical line. For smaller Hq, the constraints become more stringent. The 
upper bound on Tma x shifts towards smaller T ma x as (hfo)^- The vertical line moves towards 
smaller 7713/2 also as (hfo)^- Note that the current data prefer Hq ~ 70 Mpc/km/sec. 
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and SUSY is dynamically broken. Assuming a vanishing cosmological constant, the gravitino 
mass in this model is given by 


"^3/2 


(-^o) ^ ^DSB 

V3M, M, 


(3.1) 


In the next stage, the SUSY breaking is fed down to the messenger sector by integrating 
out the U(l)iness interaction. The messenger sector contains a gauge singlet S', whose A- and 
F-component Fs acquire VEVs after minimizing the potential.^ The scale of these VEVs are 
related to the original SUSY breaking scales Adsb as (S') (_F_g)V2 ^ 0(5(^gggADSB/167r^). 

The ratio of {Fs) to (S') determines the masses of the sparticles in the ordinary sector. By 
counting the loop factors, we obtain 

A„,ess = ^ = «^messfP\/m3/2M„ (3.2) 

where Kmess is supposed to be of 0(1) 0, and g^ess is the gauge coupling constant for the 
U(l)mess gauge interaction. In the messenger sector, there are also 5 + 5 representation of 
SU(5), i.e. SU(2)l doublets (L and L) and SU(3)c triplets {Q and Q). These have the 
superpotential 


lUniess = X 2 SLL + XsSQQ. 


(3.3) 


Once S and Fs acquire VEVs, the scalar components of L and L {Q and Q) have a mass 
matrix of the form 


/ V(F)2 X{Fs) \ 
[ X{Fs) X^S)^ ) 


(3.4) 


while the fermionic components have mass X{S). Therefore, SUSY is broken in the mass 
spectrum of vector-like 5 + 5 messenger fields. By integrating out the messenger fields, 
the soft SUSY breaking parameters are induced in the ordinary sector. With pairs of 
vector-like 5 + 5 multiplet, the gaugino masses, mci, ^ 02 , and mcs are given by 


rriGi = 


1 oF-^messYs, 

lOTT"' 


(3.5) 


where C3 = C2 = 1, and Ci = | 
sfermions (/ = mr, Jr, ^l, /'l, 


in our convention. On the other hand, the masses for the 
and cr) are given by 


m? = 2AL..A'5i:Q(Yi3) ■ (3-6) 

*One actually needs a substantially more complicated messenger sector than the original ones ill 
in order to avoid a run-away global minimum |l^. Such details are, however, beyond the scope of this 
paper. 
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Here, Ci = with Y being the usual hypercharge, and Cj = | and | if / is in the 
fundamental representation of SU(3)c and SU(2 )l, and C* = 0 for the gauge singlets. 

Combining the above relations with the experimental bounds on the sparticle masses, we 
can obtain a lower bound on A^ess, and hence the gravitino mass. A lighter state gives us 
a more stringent constraint. For most parameters, the lightest among the sfermions is the 
right handed selectron, whose mass rrie^ is given byjj] 

A^^,,A ^5 - m| sin^ 6»w cos 2f3. (3.7) 


2 f (Amess 

"" y [ 167r2 


For tan/3 close to 1, the right-handed selectron mass gives us a stringent bound on the 
messenger scale. Even if we adopt a conservative constraint of me^ > 45GeV,[| we obtain 


Amess ^ 2 X lO'^GeV X 



(3.8) 


If tan (3 is large, the D-term contribution enhances the right handed selectron mass, and 
can be larger than the experimental limit with smaller value of A^ess- In that case, however, 
the sneutrino mass 



rsj 


3 / 5'2(-^mess) 

2 V 167r2 ) 


5 / g‘f{Amess) \ 

6 y IbTT^ j 


Amess^5 + ^OS 2/3, 


(3.9) 


receives a negative contribution from the Z3-term, and Amess is still constrained to be larger 
than ~ 2 X lO^GeV. Therefore, the bound (|3.8|) holds for all values of tan/3. 

In order to translate the above constraint (|3.8|) into a lower bound on the gravitino mass, 
we need information about the gauge coupling constant of the U(l)mess; larger (/mess (at 
the messenger scale) gives us a less stringent constraint on the gravitino mass. However, if 
//mess is too large at the messenger scale, it blows up below the Planck scale or even below 
the GUT scale. I By using the 1-loop renormalization group equation, the U(l)mess gauge 
coupling constant at the messenger scale is constrained by 


tfmess 


(A„ 


< 


167r2 ~ 6mess ln(A2^t/A^g,J ’ 


(3.10) 


where A^ut is the cutoff scale where the perturbative picture may break down, and &mess = 
J2a Qmess A i^ the sum of the squared charge of the messengers. Typically, 6mess ~ 10. (For 

Tn fact, the selectron mass receives a correction from renormalization effects. However, the correction is 
less than 10%, and the following arguments are almost unchanged even if we include such effects. 

^We are aware that LEP-II has placed stronger limits on the mg, but as a function of the neutralino mass. 
The mass difference between mg and the lightest neutralino ~ i? is not large in the LEGM models, and the 
constraint weakens substantially in this situation. We chose the model-independent LEP bound because of 
this reason. 

§ There is the logical possibility of employing an asymptotically free non-abelian gauge group as the 
messenger group, and assume that its scale parameter is very close to Adsb. We, however, consider it as an 
exponential fine-tuning of parameters. 
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the model proposed in |Q, bmess = y, and for a model given in 6mess = 14.) Assuming 
bmess = 10 and Acut ~ Mqvt, 5'mess at the messenger scale is constrained to be smaller than 
~ 0.5, and 


m3/2>70keVx^ 


/ 9raess \ ^ ^ \ ^ 

V1)X/ ViSGeVy 


(3.11) 


In the minimal model, = 1, and if we assume the perturbative unihcation of the gauge 
coupling constants in the MSSM, < 4 [^]. Therefore, in any case, the lower bound 
above is about one or two order of magnitude larger than the cosmological upper bound 
(2.1). Notice that the lower bound on the gravitino mass increases as the experimental lower 
bound on the sparticle masses increases.|^ 

Based on the above estimations, we dehne the canonical set of the parameters for our 
following analysis: 


7713/2 = lOOkeV, 

(3.12) 

A^ess = 3 X lO^GeV, 

(3.13) 

{Fsy/^ = {S) = A,,ess. 

(3.14) 


Note that it is easy to raise the gravitino mass; we only have to assume a smaller value for 
the gauge coupling constant for the U(l)mess. In the following analysis, we basically assume 
the above set of parameters, and we also discuss how our results change if we vary them. 

The above estimations are based on perturbative calculations, and one may worry that 
a strong coupling in the dynamical sector may allow us to lower the gravitino mass. Such 
a scenario seems unlikely, however. To see this, it is convenient to dehne the “vacuum 
polarizations” from the DSB sector for the U(l)mess gauge multiplet: 

F.T. (0|T(A^A^)|0)^pj = iq‘^UA{q‘^)gf,u, (3.15) 

F.T. (0|T(AA)|0)ipi = (-z){#lA(g^) + ^(g^)}, (3.16) 

F.T.{0\T{DD)\0),^, = tUn{q^), (3.17) 

where F.T. stands for the four-dimensional Fourier Transform to the momentum space, and 
IPI for one-particle irreducible diagrams. At tree level, 11^ = nA = n£) = E = 0. These 
quantities are radiative corrections of if the perturbative calculation is reliable. 

The messenger scale in the LEGM model is induced by integrating out the SUSY breaking 
sector and the U(l)mess gauge multiplet. By using Ba, Ba and B/j, the SUSY breaking scalar 
mass in the messenger sector (~ the messenger scale) is given by 

~ jL../ - 4n,(,=) + nD(«=)}. (s.is) 

collorary to this analysis is that it is unlikely to have a decay of a sparticle into the gravitino inside a 
collider detector. This casts some doubts on the naturalness of ZZ77 signature at CDF in the LEGM models. 
A possible way out is to employ the vector-like model and couple a singlet field directly to the messenger 
fields in the superpotential |l^. This model, however, probably suffers from a tunneling to a color- and 
charge-breaking supersymmetric minimum if all coupling constants are 0(1) fl^. 
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If we limit ourselves to the physics of the DSB sector, there is no pole in the 11 functions at 
= 0 which, if present, implies the Higgs mechanism for the U(l)mess gauge group. The only 
singularities in H functions are, therefore, branch cuts which appear above certain threshold 
^ Apgg which is the only scale in the problem. Then the integrations in Eq. ( p.l8| ) can 
be Wick rotated and we obtain 




jmess 

Ibvr^ 


roo 

/ dql{3IlAi-ql) - AIlx{-ql) + Uoi-ql)}. 

JO 


(3.19) 


Now it is clear that the integration is purely Euclidean, and hence all H functions are far 
off-shell. Thus, the perturbative result is essentially reliable even when the DSB sector is 
strongly coupled. It is also useful to recall that similar calculations of vacuum polarization 
amplitudes in QCD tend to agree with lowest order perturbative results for the running of 
the hne-structure constant, or the scaled-up QCD estimate of the electroweak S'-parameter.|U] 
We therefore conclude that there is no signihcant enhancement of the resulting and 

hence the estimates of the messenger scale and the resulting gravitino mass ( p.ll|) can be 
trusted. 

The constraint (|3.11|) sets severe bounds on cosmology. In particular, we need some 
mechanism to generate a dilution factor of ~ (m 3 / 2 / 2 h^keV) at a relatively low temperature 
below the upper bound on the maximum temperature given in Eq. (|2.^) and in Fig. 2, if the 


gravitino mass is larger than 2h^keV. Furthermore, even if we adopt such a large entropy 
production at a low temperature, baryogenesis may still be a problem. The Affleck-Dine 
mechanism [jTI[ for baryogenesis is one of the possibilities to generate baryons at a relatively 
low temperature.^ 

However, in the LEGM models, the behavior of the flat direction at large amplitude is 
quite different from the usual supergravity case. Thus, even if we assume the Affleck-Dine 
mechanism, it is a non-trivial question whether we can have enough baryon number density. 
In the following sections, we pursue this possibility, and as a result, we will see that the 
Affleck-Dine mechanism works sufficiently well, enough to explain the present value of the 
baryon-to-entropy ratio. 


II Note that the S'-parameter is defined by the vacuum polar ization amplitudes at = 0, and hence more 
sensitive to the non-perturbative effects than our case ( |3.19 ) which smears them over the wide range of 
q^. Still, a perturbative estimate of S differs from the scaled-up QCD only by a factor of two. One may 
also estimate the S'-parameter by assuming that it is dominated by the p and ai poles. Then the result is 
obtained by the tree-level process. Even so, the coupling of the resonances to the current operator has a 
factor of l/47r and the counting of l/47r factors remains the same in as the perturbative one-loop result. HM 
thanks Christopher D. Carone on this point. 

** Electroweak baryogenesis may be another possibility to generate baryon asymmetry at relatively low 
temperatures. However, the resulting baryon-to-photon ratio depends on the details of the complicated 
dynamics of the phase transition. Furthermore, the generated baryon asymmetry would not be large, if any, 
and would probably not survive a huge entropy production to dilute the string moduli fields as discussed in 
Section 
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4 Flat Directions in the LEGM Models 


As discussed in the previous sections, the constraint from the gravitino cosmology is quite 
severe in models with the LEGM. Therefore, it is preferable to look for baryogenesis scenarios 
which do not require high temperatures. 

We will focus on the Affleck-Dine baryogenesis in the LEGM models in this paper. There 
is one crucial difference from the hidden sector models: the potential along the MSSM flat 
directions is not simply parabolic. Therefore we discuss the form of the potential hrst in this 
section, and then its implication to the Affleck-Dine baryogenesis in the next section. 

In the hidden sector models where the supersymmetry breaking effect is mediated by 
Planck-scale operators, the soft supersymmetry breaking parameters are actually “hard”, in 
the sense that they renormalize as usual mass terms between the Planck-scale and the weak 
scale. On the other hand, the supersymmetry breaking scalar mass terms are suppressed 
beyond the messenger scale in the LEGM models. 

This is analogous to the situation in the QGD. The current masses of the quarks renor¬ 
malize according to the ordinary perturbation theory. They are “hard” masses. However the 
constituent quark masses are suppressed by a power of the energy: “soft”. This is because 
the constituent quark masses are dynamically generated by the spontaneous chiral symmetry 
breaking, which is characterized by the order parameter {qq). The constituent quark mass 
has to be proportional to this order parameter. At high momentum transfer, a dimensional 
analysis then tells us that the effective constituent mass behaves as mconst(Q^) 

The same argument applies to the soft supersymmetry breaking masses from the LEGM. 
Supersymmetry is broken by an F-component of a chiral superheld, (F 5 ) 7 ^ 0. The soft 
supersymmetry breaking scalar mass is necessarily proportional to the order parameter of 
supersymmetry breaking, he., rn? oc {Fs)\Fs). A dimensional analysis tells us that it is 
suppressed at high momentum transfers, m?{Q‘^) ~ {F{Fs)/. Therefore, the supersym¬ 
metry breaking parameters “shut off” at high energies. 

The potential of a MSSM flat direction is given simply hy V = where is a 

soft supersymmetry breaking mass. A renormalization group improvement gives us H = 
m^(|0p)|0p. In the hidden sector case, m^(|0p) has only a logarithmic dependence on |0p 
and hence can be taken approximately constant unless it crosses zero at some energy scale. 
For most cosmological applications, this is a sufficiently good description. In the LEGM 
models, however, the effective mass exhibits a power dependence on |(/)p which 

cannot be neglected. We expect that behaves as {Fs^(Fs) for large |0|, and 

hence the potential behaves approximately like a constant for \(f)\ > {S), which is the mass 
scale of the messengers. 

We performed an explicit two-loop calculation of the effective potential V (0) and its 
details are presented in Appendix Here we only quote the result. As expected, the 
potential behaves parabolically around the origin, while it becomes approximately constant 
for large |0|; actually it keeps growing slowly as (ln|0p)^. The potential of a MSSM flat 
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direction behaves as 

for small |0| -C (S), and 

(4.2) 

with 

^0 ~ ~ (3 X 103GeV)^ (4.3) 

for large |0| 3> (S'). Here, g generically refers to standard model gauge coupling constants.^] 
For extremely large 101, however, the contribution from supergravity becomes important. 
Supergravity generates a contribution to the scalar potentials ~ m3y2l0P for any |0|.[| To 
determine the relative importance of the LEGM and supergravity contributions, we compare 
their derivative V because this is the quantity which appears in the equation of motion. The 
derivative of the potential from the LEGM is 


dV 20* 

84, °\4\^ 



(4.4) 


which is to be compared with the supergravity contribution dV/dcj) = The super¬ 

gravity contribution is more important above a threshold value 0eq which is given by 


0eq ~ 


2Ho 

^ 1/2 



1/2 

~ 7 X 10^^ GeV X 


/ ^ 3/2 

UOOkeVy 




1/4 


3 X 103 GeV 


(4.5) 


The motion of the flat direction is determined by the effective potential given in this section 
and the canonical kinetic term, and there is no need to include the wave function renormal¬ 
ization factor at this order in perturbation theory. See Appendix ^ for details. 


5 Affleck—Dine Baryogenesis in the LEGM Models 


In this section, we aim at estimating the size of baryon-to-entropy ratio from Affleck-Dine 
baryogenesis in the LEGM models. Because of the multiple scales in the problem, the 


* The readers may wonder why Eqs. (^) and (O) have only two powers of gauge coupling constants 
despite the two-loop-ness of the effective potential. This is the result of an explicit calculation, and we can 
also explain it in a simple way. When the field value is large, the standard model gauge multiplets acquire 
large masses of order g(j>. The effective potential is generated by the exchange of heavy gauge multiplets, 
and hence it is suppressed by l/|g0p. This cancels two powers in gauge coupling constants. 

Mhis is true for the minimal supergravity and its variants. This contribution, however, is not there in 
no-scale supergravity, or in general, supergravity with Heisenberg symmetry In such a case, 0eq 

must be taken at M* in the rest of the discussions. 
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discussion becomes somewhat complicated. The basic conclusion is that the Affleck-Dine 
baryogenesis works as efficiently as in the hidden sector case, but in a much more non-trivial 
manner. Finally we discuss the possible dilution of gravitinos via the decay of the Affieck- 
Dine flat direction, and hnd that the gravitinos can be diluted below the closure limit if the 
initial amplitude of the flat direction is sufficiently large. 


5.1 Generalities 


In Affieck-Dine baryogenesis [j^, one assumes that a MSSM flat direction has a large am¬ 
plitude at the end of the primordial inflation. The mechanism to achieve a large amplitude 
varies: a negative curvature from a non-minimal Kahler potential |^, or no-scale super¬ 
gravity [^. In any case, it tends to be equal to or larger than the expansion rate of the 
Universe during inflation ilinf ~ 10^^ GeV - 10^^ GeV depending on inflationary scenarios. 
We phenomenologically parameterize it just as the initial amplitude 0o- 

A typical assumption is that there is a baryon-number violating Kahler potential term 
K ~ rq*u'^d^/Ml, where M* is the reduced Planck scale.|3 The supersymmetry breaking 
effects from the LEGM generates a term in the potential^ 


O 



9 


\e‘^Fs\^ 


(47r)4 

2 \ 2 


In 


(^)V Ml 


rq*u^F + h.c. 




(5.1) 


while the supergravity effect induces an operator in the scalar potential 






Ml Ml 


l*q*u^(f + h.c. = -Mj^rq*u'^(f -f h.c. 


Ml 


(5.2) 


The LEGM operator is dominant if \(j)o\(j)eq, while the supergravity one dominates if 
|0o|~0eq- Therefore, we discuss the two cases separately below. In either case, the size 
of the baryon-number violating operator is much smaller than in the hidden sector case 
(Appendix 0. It turns out, however, that the baryogenesis proceeds efficiently with these 
operators. 

Below we generically refer to the helds as 0 without distinguishing various species. The 
baryon number in the scalar sector is given by 


ub = (5.3) 

*It could also well be the GUT-scale Mgut- However, in this paper, we assume this form of the baryon 
number violating operator for simplicity. The extension to the case with K ^ is trivial, and 

one can easily estimate the resulting baryon-to-entropy ratio. 

^We have not calculated this explicitly. This form is expected based on the analogy to the calculation of 
the effective potential in the previous section. The only difference is that the previous one arises from the 
kinetic term (/)*d in the Kahler potential rather than from a non-renormalizable term l*q*u'^F/Ml here. 
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while the baryon number violating operator is written as 


O ~ 



1 


(0" + n- 


(5.4) 


5.2 |0o| 


> 

rsj 


^eq 


For sufficiently large |0o| (|0o|~0eq), the supergravity contribution is initially important, 
and the field begins to roll down the potential when the expansion rate of the Universe H 
is comparable to Ff ~ '^ 7 . 3 / 2 . Let us first estimate the primordial baryon number asymmetry 
in this case, |0o| ~ 0eq- 

A rough estimation of the baryon asymmetry, which is generated just after the start of 
the oscillation of the 0 field can be done only by using simple order of magnitude arguments. 
With the above baryon number violating operator, the time evolution of the baryon number 
is given by 

Ab + mn, = i ■ (5.5) 

When the field begins to roll down the potential, its initial motion is slow, and one can 
neglect the ub term in the equation (see Appendix 0). Then the resulting baryon number 
can be estimated by |2^ 

t (dO^ dO \ 

~ ■ ayU ■ 

Using the approximate order of magnitude of the operator and H ~ 7713 / 2 , we obtain 


Ub 


7?73/2lm(0g) 

M2 


(5.7) 


It depends on the imaginary part of the initial amplitude. The entropy of the radiation at 
this stage is given by s ~ while the energy density prad ~ g*T^ ~ ^ 3 / 2 ^*- ^y putting 
them together, we can estimate the baryon-to-entropy ratio. 


Ub 

s 


rsj 



lm(0^) 


~ 4 X 10^° X 



f_rri3/2_\ 

UOOkeVy 


-1/2 

sin 46*0, 


(5.8) 


where the initial amplitude is parameterized as 0o = As one can see, a large baryon 

asymmetry can be generated, if the initial amplitude of 0 is not too small. Therefore, a large 
enough baryon number can remain in this scenario, even if there is a substantial entropy 
production after the Affieck-Dine baryogenesis. Note that one obtains exactly the same 
expression in the hidden sector models, but with a different 7773 / 2 . 

The present baryon-to-entropy ratio is also given by the above formula, if there is no 
significant entropy production. However, in a realistic situation, there can be entropy pro¬ 
duction. In particular, the decay of the Affieck-Dine field 0 may produce a large amount 
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of entropy. Furthermore, for m 3/2 ~ keV, we need a non-negligible entropy production 
to dilute the primordial gravitino. If there is an entropy production after the Affleck-Dine 
baryogenesis, the primordial baryon number density is also diluted. In the following, we 
discuss how an entropy production affects the results. 

We estimate the entropy production due to the decay of the flat d irection. As discussed 
in the previous section, 0 starts to oscillate when T = Tq ~ if |0o| ^ 0eq- 

During |0|^0eq5 the potential for 0 is dominated by the supergravity contribution, and 
hence 


= const. <l<#'l<l'#’ol). 


(5.9) 


Thus, the temperature at 


0eq, which we denote Ten, is estimated as 


2/3 




Tn 


^eq 


\‘P« 


(5.10) 


For 101 ^ 0eqt the potential for the flat direction is dominated by the LEGM piece, and the 
evolution for 0 does not obey the relation ( |5.1j| ). By using the virial theorem, the evolution 
of the flat direction can be estimated, and is given by 


|0|T^ = const. ((S') < 




^eq 


(5.11) 


(See Appendix 0.) By using the above relations, we obtain the dilution factor due to the 
decay of the flat direction. 

Now an important question is at what held amplitude 0 decays into radiation. When the 
motion is dominated by a parabolic term, we know the time dependence of the oscillation 
(just a harmonic oscillator), and we can calculate the rate of particle production in such 
a background. The result is known to be the same as the single particle decay rate, if 
the amplitude is not too large compared to the oscillation frequency. Once the amplitude 
becomes comparable to (S'), the potential is almost parabolic, and we hnd that the coherent 
oscillation decays into radiation rapidly. On the other hand, a corresponding calculation 
is difficult when the potential is dominated by the logarithmic term. In our analysis, we 
regard the decay amplitude, 0 dec 5 as a free parameter, and discuss the 0dec-dependence of 
the results.fi 

The temperature of the background radiation at the decay time of the flat direction, Tdec, 
is given by 


Tdec 




(5.12) 


^For a canonical parameter, we take 0dec 10® GeV for estimating dilution factors. Since 0 decays at 
0 ^ {S) ^ 3 X 10"^ GeV at latest, this choice gives us the minimum estimate of the baryon asymmetry. If 
0 decays earlier, the dilution factor is less and the baryon asymmetry is larger. We likely overestimate the 
dilution factor with this choice. We will come back to this point later when we discuss a possible dilution of 
gravitinos from the decay of the flat direction. 
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On the contrary, the energy density of the flat direction is 


Pflat ~ ho- 


(5.13) 


Then, if pfjat ~ Prad, the dilution factor from the decay is given by 


D 


\Prad/ \fi'*^dec/ 


3/4 


6 X 10® 


a/4 


I0o| \ f 0dec \ / ai3/2 \ ^ 

mJ l^lO^GeV; UOOkeV; IsxlO^GeV 


(5.14) 


Note that H ~ 1 if poat ~ Prad, or in terms of the initial amplitude, 0o ~ 10^“^ GeV. 

Gombining the above dilution factor with the estimation of the primordial baryon number 
density given in Eq. (|5.8|) , we obtain the present baryon number asymmetry. In order to 
make a pessimistic estimate of the resulting baryon asymmetry, we assume that the flat 
direction decays only when its amplitude is as small as (S). This assumption maximizes the 
entropy production, and hence, gives us the minimum value for the baryon asymmetry. If it 
decays earlier, then the entropy production is less and hence the baryon asymmetry is larger. 
With this caveat in mind, we can make an estimate of the resulting baryon-to-entropy ratio, 
and in the case with entropy production {D > 1), the resulting baryon-to-entropy ratio is 
given by 


Ub 

s 


1/4 




~ 70 X 



/ ^dec ^ 

1^105 GeVy 




1/4 


3 X 103 GeV 


sin 400. 


(5.15) 


Note that the result is independent of the gravitino mass. It is intriguing that the final result 
is more or less the same as in the hidden sector case Eq. ( P7.9|) . 


5.3 1001 ~ 0eq 

Next, we discuss the case of |0o| ^ 0eq-0 fhis case, the potential for the flat direction is dom¬ 
inated by the LEGM-piece, and the flat direction starts to move when H Vl^'(</>o)|/|0o| ~ 
Y^Vo(ln |0oP/(5')^)/1001- The temperature at this stage. To, is estimated as 

§As noted before, the no-scale supergravity does not generate potential term proportional to 1713^2 
hence the evolution of the flat direction is always dominated by the LEGM piece. Then the formulae 
presented in this subsection must be used even for a larger |(/io| ~ Such a large |^| is indeed expected 
in the no-scale case because the flat directions remain flat even during the inflation ||^. 
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Then, by using Eg. (^-61 ), we can estimate the resulting baryon number density, and hence 
the baryon-to-entropy ratio. Note that the baryon-number-violating operator O is different 
from the previous case. 

Following exactly the same steps as in the previous case, we find 


Ub ~ 





3/2 


V«'^\4'o? 

Ml 



3/2 

sin 46*0, 


(5.17) 


when the flat direction starts to move. Therefore, the baryon-to-entropy ratio is given by 


Ub 

s 


7-1/4. 




In 




o| 


12\ 3/4 




sin 46*0. 


~ 6 X 10 "^^ X 


100 


9/2 


En 


1/4 


1012 GeV 


3 X 103 GeV 


sin 46*0 


(5.18) 


which is typically too small. It is useful to note that the above formula is larger by a factor 
of the logarithm than the corresponding formula Eg. ( b.8|) for the case 0o ~ 0eq when we 
substitute 0o = 0eq- Of course such a discontinuity cannot exist. It simply means that there 
is a transition region at 0o ~ 0eq where there is a slight rise in when we cross 0o ~ 0eq 
from above. 

The dilution factor can be estimated also along the lines of the previous case. We have 


^dec ~ To 



(5.19) 


and hence (if peat > Prad), 


D ~ 




Pradj ~\a.TL:) ^ 


2 X 10"^ 


( I0ol \ 

1^1012 GeV; 


5/2 


/ 0 dec \ 

1^105 GeVy 


(5.20) 


(More correctly, the dilution factor is H = ((pfiat + Prad)/Prad)^'^"^ and cannot be less than 
unity.) Therefore, the dilution factor is much less important than in the previous case. 


5.4 Numerical Analysis 

A more detailed behavior of the baryon-to-entropy ratio can be studied by a numerical 
calculation. Here, to see the behavior of the results, we show the resulting baryon-to-entropy 
ratio for a particular set of parameters: Vq = (3 x lO^GeV)^, (S) = 3 x 10“^ GeV, 0dec = 10® 
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GeV, and sin46*o = 1. The results for other sets of parameters can be easily estimated by 
using Eqs. (l^) and (|5.15|) . First, we solved the equation of motion for the flat direction, 


• dV 

4> + 3H4> + — = 0, 


(5.21) 


with the potential^] 


V = m 


3/2 I 


+ Go In 




Gn 


2 I ’'0 

7713/2 + 


In 


i0in (5) 


M2 


( 0 " + 0 ^ 


*4\ 


(5.22) 


We start our calculation at a temperature much higher than Tq (T = IOTq), and follow the 
evolution of the flat direction as well as the temperature of the thermal bath. With the 
initial value 0o = 0.2M*e*’^/®, the initial motion is shown in Fig. As one can see, 0 starts 
an elliptical motion due to the baryon-number violating term in the potential. This means 
that a non-vanishing baryon number is generated once 0 starts to oscillate. We indeed found 
the generated baryon number to be consistent with our estimates in the previous subsections 
within a factor of a few. 

With the motion of the flat direction, we can calculate the baryon number by using 
Eq.(|OD, and hence the baryon-to-entropy ratio. After several cycles of oscillation, the 
baryon-to-entropy ratio becomes almost constant. Then, the evolution of 0 and T can 
be easily traced by using Fqs. (|5.9|) and ( |5.11|) with RT = const. Finally, we calculate the 
dilution factor, D = ((prad + P</.)/Prad)^^^, at the decay time of 0, and multiply the primordial 
baryon-to-entropy ratio by D~^ to obtain the resulting baryon asymmetry. 

In Fig. we show the |0o| dependence of the present baryon-to-entropy ratio, ub/s. 
From the hgure, we can see that the results based on the order of magnitude estimations 
provide us good approximations. For a sufficiently large 0o such that the entropy production 
is signihcant, the resultant baryon-to-entropy ratio is independent of the gravitino mass, and 
is proportional to |0oP- We also checked that the approximate formula (|5.15|) reproduces 
the behavior for the large |0o| region. For a smaller value of |0o| ~ 10^^ GeV, the entropy 
production from the decay becomes negligible. The result then is proportional to |0o|^ 
(Fq. (^.8|)). For an even smaller |0o| ^ 0eq ~ 7 x 10^^ GeV x (100 keV/7?73/2), the behavior 
goes over to |0o|^^^ (Eq. (|5.18|) ). As noted in the paragraph below Fq. (|5.18|) , there is a 
transition region from ~ lO0eq to 0eq, where the curves fall less steeply because a logarithmic 
enhancement factor comes in. In any case, the baryon-to-entropy can clearly be sufficiently 
large in this scenario, as required by the standard big-bang nucleosynthesis rr^/s ~ 10“^°, 
if the initial amplitude is larger than 10^^“^^ GeV. 

A more precise estimate of the baryon asymmetry requires the specihcation of the flat 
direction, the relevant operator, and the size of the initial amplitude. The usual caveat 
concerning the B — L invariance applies: If one employs an operator which preserves B — L 


^The potential given in Eq.( b.22|) is unbounded-below for and higher-dimension operators are 

supposed to stabilize it. However, we only consider the initial amplitude to be less than ^ M* in our analysis, 
and hence the postulated potential is a good enough approximation. 
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Figure 2: The initial motion of the flat direction with the potential given in Eq. 
Here, we take 7713/2 = 100 keV, and 0o = 0.2M*e*’^/®. 
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Figure 3: The resulting baryon-to-entropy ratio as a function of the initial amplitude 00- 
The parameters are taken to be Vq = (3 x lO^GeV)^, Oq = vr/S, ^dec = 10® GeV, and 1713/2 = 1 
keV (dotted line), 1713/2 = 100 keV (solid line), and 1713/2 = 10 MeV (dashed line). 
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symmetry like that generated by SU(5) grand unified models, it may be wiped out again by 
the electroweak sphaleron effect . One needs to either preserve the B and L asymmetries 
using a Bose condensate |^, or by generating a. B — L asymmetry |]^ , |^ . In view of the 


discussions of the next section, we find the protection of B via a Bose condensate is a likely 


scenario. 


5.5 Diluting Gravitinos 


In the LEGM models, the mass of the gravitino is about 100 keV, and its mass density 
exceeds the closure limit if T^ax is larger than (1 - 10) TeV, as discussed in Section ^ 
Since the reheating after the primordial inflation raises the temperature typically above 
^RH ~ 10® GeV or so, we assume that the gravitinos were once thermalized, which is of 
course the worse case scenario. We will discuss whether the decay of the flat direction can 
generate a large enough entropy to dilute gravitinos below the closure limit. 

Before discussing the implication of the entropy production to the gravitino, it is useful 
to estimate the freeze-out temperature of the gravitino, T^-^eezei which is the temperature 
at which the expansion rate of the Universe becomes comprable to the production rate of 
the gravitino: if(Tfi.eeze) ~ Stot’^rad(^freeze)- By usiug the gravitino production cross section 
given in Eq. (|2.5|) , we obtain 


^freeze ~ 200TeV X 


/ m^/2 Y 
UOOkeVy 


/ mG3 \ ^ 

UTeVy 


(5.23) 


Below Tfreeze; the expausiou rate of the Universe becomes larger than the production rate 
of the gravitino, and hence the gravitino cannot be thermalized. Therefore, an entropy 
production at T ^ Tfreeze dilutes gravitinos produced before the entropy production. On the 
contrary, even if the entropy is produced when T ^ Tfreeze, the gravitino is thermalized again, 
and its number density is determined by the thermal distributions. If the energy density of 
the flat direction dominates the energy density of the Universe, decay of the flat direction 0 
reheats the Universe. The reheating temperature is estimated as Tr (Uo)h4^ 1-10 TeV. 
Gomparing this reheating temperature to Eq. (|5.23|) , we can see that the gravitino cannot 
be thermalized after the decay of 0. In other words, gravitinos produced before the decay 
of (j) are diluted with a dilution factor given in Eq. (|5.14|) , if the decay of (j) produces the 
entropy. 

If the gravitino mass is larger than ~ 1 keV, we need a substantial entropy production; 
otherwise, the Universe is overdosed by the mass density of the gravitino. By assuming that 
the gravitino is thermalized, we can estimate the number density of the gravitino as 


^ 3/2 (T) 


3 g*{T) 

2 (7* (^freeze) 


Wad(^)- 


Then, requiring fl ^/2 = D ^m 3 / 2 %/ 2 /Pc < 1, we obtain 


d > 50 X h-2 


/ m3/2 \ 

UOOkeVy 


(5.24) 


(5.25) 


20 


















Comparing the above constraint with Eq. ( |5.14|) , we can see that the decay of the flat 
direction can produce enough entropy to dilute the gravitinos away. For example, for 0dec ~ 
10^ GeV and m 3/2 ~ 100 keV, the dilution factor is large enough, if |0o| ~ 10^"^“^® GeV. Even 
with such a large dilution, the present baryon-to-entropy ratio can be sufficiently large (see 
Eq. (|5G5|) ). 

If the flat direction decays at an amplitude larger than ~ 10^ GeV, the dilution factor 
given in Eq. (|5.14|) becomes smaller, and the entropy production due to the decay of 0 
may not be enough to decrease the gravitino density. In addition, for 7 / 13/2 ~ 100 keV, 
the reheating temperature Tr 1 — 10 TeV may be higher than the freeze-out 

temperature of the gravitino. In these cases, we have to assume an extra source of the 
entropy production of O ( 7713 / 2 /2/i^keV) after the reheating. Even in this case, the estimation 
of the primordial baryon-to-entropy ratio (|5.8| ) is still valid, and the hnal baryon asymmetry 
can be as large as the one estimated by Eq. (^.15|) and an additional dilution factor Eq. (|5.25|) 
required to dilute the gravitinos. Therefore, Affieck-Dine baryogenesis can generate a big 
enough baryon asymmetry to explain the present value of the baryon-to-entropy ratio. 

In fact, we can crudely estimate the decay amplitude 0dec even when it is in the region 
of the logarithmic potential, if its motion is circular rather than elliptic. The change from 
the original Affieck-Dine estimate of the decay rate Eq. ([E.7|) is that the rotation frequency 
of the (j) held is given by (Vq ln(|0p/(S')^))^/^/|0| rather than 7713 / 2 . Since we are interested 
in a dilution factor, we assume that the hat direction dominates the Universe, and the held 
decays when 


— 


a. 


Voln(|0dec|V(^)')\'^' 


' 1 4’dec I ^ 


|0decP 


H 


r\j n r\j 


|0d '' 


and 

|0dec| ~ 




2\ 1/2' 


1/5 


M* 


In 


2dec I 




(5.26) 


8 X 10^ GeV X 


Un 


(3 X 103 GeV) 



Therefore, the decay amplitude does not change much from the value assumed before. On 
the other hand, the case with an elliptic orbit is more difficult to deal with. We are not aware 
of any study on the decay rate of 4 for an arbitrary elliptic motion even for the parabolic 
potential. The other limit of almost linear motion is discussed in the literature and tends to 
give a larger decay rate, and hence a larger 4 de 


28 


However, we believe the motion of the 
4 held in the case of our interest here, namely for |0o| ~ GeV to dilute gravitinos, to 

be quite far from a linear one, and 4 dec is not likely to be much larger than our estimate. 
We conclude that 4dec is not much larger than the minimum possible value lO® GeV which 
we used in most of our discussions. 
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6 Cosmology of String Moduli 


We point out in this section that the moduli helds in the string theory, if they acquire masses 
in the LEGM models, are stable and drastically overdose the Universe. 

According to a general analysis [^, string moduli acquire masses comparable to the 
gravitino mass m 3 / 2 . Their initial amplitude is likely to be of the order of the string or Planck 
scale because it is the only scale in the problem. The cosmological problem of the moduli 
helds is discussed extensively in the literature in the context of hidden sector models (for the 
original paper, see |^). There, the moduli helds acquire masses of the order of 1 TeV, and 


decay after nucleosynthesis, thereby spoiling the success of the nucleosynthesis theory. Even 
if one pushes the mass to 10 TeV so that the moduli helds decay before nucleosynthesis, the 
enormous production of entropy with a dilution factor of order M*/m 3/2 ~ 10^^ wipes out all 
pre-existing baryon asymmetry. This problem may be solved by adopting the Afheck-Dine 
baryogenesis |^, or by the thermal inhation |]^ (see Appendix ^ for more discussions). 

In the LEGM models, the situation is completely diherent.|^ The string moduli are stable 
within the cosmological time scale, and are still oscillating around their potential minima. 
A dimensional analysis gives the decay rate of a moduli held to be T ~ ~ 

(3 X 10^® years)“^, for 777 . 3/2 ~ 100 keV. Therefore there is a problem concerning its energy 
density. 

The estimation of the moduli energy density is straight-forward. When a moduli held 
begins to oscillate, the expansion rate is TT ~ 7773 / 2 . The entropy at this stage is given by 
s ~ (7y^(7773/2M*)®/^. Assuming the initial amplitude of order M*, the ratio of the moduli 
energy density to the entropy is given by 


Pmoduli 


~ 1.3 x 10® GeV 


»^ 3/2 \ 

100 keV/ 


( 6 . 1 ) 


Since both the energy density of the moduli and the entropy are diluted by the expansion with 
the same rate the ratio remains constant until now unless there is entropy production. 
On the other hand, the total energy density is bounded from above by the critical density 

Pci 

< _^ = 3.6 X 10 - 9/72 GeV. (6.2) 

^ '^now 

where 5now is the present value of the entropy density. The predicted ratio Eq. o is in 
gross conhict with the constraint Eq. (|6.2|) . 

It is not clear how such an enormous energy density can be diluted. First of all, the 
necessary dilution factor is at least 10^®. Furthermore, one needs such an entropy production 
at a very late stage of the Universe, with H < 7773/2, or equivalently, T < ( 7773 / 2 = 

1.3 X 10® GeV for 7773/2 = 100 keV. One needs to create a baryon asymmetry either after 

*It was argued that the problem does not exist if SUSY is broken dynamically, which is true for scalar 
fields which directly participate in the dynamical supersymmetry breaking. However, the string moduli fields 
were not considered in this discussion. 
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such an enormous entropy production at a very late stage, or large enough to survive the 
enormous entropy production. 

Actually, the Affleck-Dine mechanism may create a large enough baryon asymmetry to 
survive the enormous entropy production which dilutes the string moduli below the critical 
density as we will show below. 

A quantity which remains constant over an entropy production is the ratio of the baryon 
number and the moduli energy density, because both of them scale as R~^. We estimated 
in the previous section that the initial baryon number density is ub mg/alm^^/M^, at the 
time when the flat direction begins to oscillate, he., H ~ 7713 / 2 - In fact, this is the same 
time as when the moduli helds begin to oscillate, and the energy density of the moduli is 
Pmoduii ~ '^ 3 / 2 ^* - Therefore, their ratio is determined at this stage: 


Pmoduli 

- ~ 1713/2 

Ub 



(6.3) 


or equivalently. 


Pmoduli _i /|0o|\ . ^/i 

X"^3/2(^J sm40o. 


s s 

Combining the above equation with the constraint 

4 

4 X 


^B ^ Pc -1 

— <-X m 


\(t>, 
l M. 


0 


sin 401 


0 


I, we hnd 
"13/2 


100 keV 


(6.4) 


-1 


100 




sin40n. 


(6.5) 


As one can see, if |0o| ~ 10^^ GeV, the baryon-to-entropy ratio may be larger than ~ 10“^° 
even if we assume a large entropy production to dilute the moduli held. 

The important question is whether one can have a brief period of inhation at such a late 
stage of Universe to dilute string moduli in the LEGM models. The inhationary expansion 

100 keV, he. the energy density of the inhation 


rate must be less than 7713/2 


Pinf ~(m3/2M*)2 (10' GeV)^. Moreover, the e-folding should not exceed 20 or so in order to 

keep the primordial density huctuation generated by a “standard” inhation with if inf ~ 10 


11 


10^^ GeV [^|. On the other hand, an e-folding of iV ^ 5 is sufficient to dilute the string 
moduli by 10“^^. A thermal inhation may oher a natural solution to these questions. 

Fortunately, we apparently do not need to introduce new energy scales into the model in 
the framework of the thermal inhation. Suppose a positive mass squared of m? ~ (100 GeV)^ 
is generated for a scalar held y due to higher order loops at the energy scale of Adsb ~ 
10^ GeV. The renormalization group running of the mass squared may drive it negative at 
a scale slightly below A^gB. If the scalar held is a hat direction of both F- and H-terms in 
the potential, it develops a minimum at Absb. This is an ideal potential for a thermal 
inhation. The scalar held may initially be stuck at the origin because of the thermal ehects, 
giving a cosmological constant. As the radiation gets red-shifted, the thermal ehects turn 
oh and the held rolls down the potential to its true minimum y = u. The e-folding in this 
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case is roughly ~ ^ ln(u/m) In(ADSBM) 5 []^ which is exactly what is needed to 
dilute the string moduli below the critical density. 

It may be interesting to compare this result with the case of the hidden sector SUSY 
breaking scenario with the Polonyi field or with the string moduli. (Hereafter, we call them 
generically as “Polonyi fields”. See Appendix ^ for the estimation of baryon asymmetry in 
this case.) Note that Eq. (|6.4|) is valid. Thus, the question is the constraint on the energy 
density of the Polonyi held, pz- In this case, the Polonyi held decays much faster than the 
LEGM case since its mass is larger, and a typical lifetime for the Polonyi held is given by 


'N, 


Tz 


ch "^3/2 


dvr M? 


-1 


lO^sec X 

VlTeVy ’ 


( 6 . 6 ) 


where Nch ~ 0(10) is the number of the decay channel. Thus, it does not contribute to the 
mass density of the present Universe, and the constraint (|6.2|) cannot be applied. However, 
it may ahect the great success of the standard big-bang nucleosynthesis (BBN) scenario. 

The mass density of the Polonyi held speeds up the expansion rate of the Universe when 
the neutron decouples from the thermal bath (he., T ~ IMeV), which may result in an over 
production of ^He. Furthermore, the radiative decay of the Polonyi held induces cascade 
photons which cause the photohssion process and change the primordial abundances of the 
light nuclei. The constraint on the primordial density of the Polonyi held strongly depends 


on its lifetime [^. If ^ 10^ sec, nucleosynthesis requires Pz/s^ 10 ® GeV. For a Polonyi 


held with a longer lifetime, the constraint becomes more stringent. For a Polonyi held with 
10 ( 4 - 5 ) 

sec, which is the case for the Polonyi mass typically larger than a few TeV, its 




< 


mass density is constrained as Pz/s^ 10“^ GeV. These constraints on p^/s are compared to 
the estimate of the baryon-to-entropy ratio (|E.14|) which holds irrespective of the presence 
of a substantial dilution of Polonyi held by, e.g.^ a late inhation. 


ns Pz -1 



sin 400. 


Thus, for this range of the Polonyi mass, the resulting baryon-to-entropy ratio may still 
be as large as lO^^o if |0 q| ~ M*, and hence the Afheck-Dine scenario may provide us a 
reasonable value for the baryon-to-entropy ratio. However, if the Polonyi held has a longer 
lifetime, as for a sub-TeV Polonyi mass as usually expected, the constraint on pz becomes 
even more stringent. In particular, for the case with Tz ^ 10^ sec, which typically corresponds 
to ruz^ 100 GeV, Pz/s^ 10“^3 GeV. In this case, the result is too small to be identihed as 
the present baryon asymmetry of the Universe. 


7 Conclusion 

We studied the cosmology of the LEGM models. We hrst estimated the lower bound on 
the gravitino mass, and saw that the bound conhicts with the cosmological constraint if 
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the primordial gravitino is not diluted. This fact indicates a huge entropy production at 
a relatively low temperature, and the conventional scenario of baryogenesis may not work 

well. 

In this case, the Affleck-Dine baryogenesis is one interesting possibility. The size of the 
baryon number violating operators is much smaller than in the hidden sector models. How¬ 
ever the flat direction begins to move at a much later stage which in turn increases the baryon 
number. The dilution factor due to the decay of flat direction also has a complicated depen¬ 
dence on parameters. After putting all the effects together, we found that the Affleck-Dine 
baryogenesis works efficiently for an initial amplitude of the flat direction, |0o| ~ 10^^ GeV. 
We also discussed that the decay of the MSSM flat direction may provide enough entropy 
to dilute the primordial gravitino for a relatively large initial amplitude of the flat direction, 
|0o| ~ GeV. Therefore, the gravitino problem in the LEGM models may be solved if 

we assume such a large initial amplitude. 

We also discussed the cosmological implication of the moduli fields in the string theory. 
Their masses are of the order of the gravitino mass, and their lifetime is much larger than 
the present age of the Universe in the LEGM models. The mass density of the moduli held 
may overdose the Universe. To dilute the moduli helds, a very late inhation is needed. 
We found that the baryon asymmetry generated by Affleck-Dine baryogenesis can be large 
enough to survive such a late inhation for |0o| ~ 10^^ GeV, even if we assume a huge entropy 
production to dilute the primordial moduli held below the critical density. 
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A Computing the Effective Potential 

The ehective potential for the hat direction 0 can be computed by the following usual 
procedure. We let it to have an expectation value, and calculate the vacuum energy in the 
presence of 0 background. The vacuum energy is identically zero if we do not pick the ehect 
of supersymmetry breaking. The lowest order contribution is from two-loop diagrams where 
the standard model gauge multiplets couple to the vector-like messenger helds whose mass 
spectrum breaks supersymmetry. The gauge multiplets acquire masses because of the 0 
background and hence the result depends on 0. 

The mass spectrum of the messenger sector is M for fermions, and ± MB 

for scalars. We assume one vector-like multiplet with a unit U(l) charge, and calculate 
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Figure 4: Feynman diagrams which contribute to the vacuum energy in the background of 
the flat direction 0 = 0. The vertices are due to the D-term potential. The scalar held with 
mass m = 2g{(j)) is the scalar component of the massive gauge multiplet in the presence of 
the background 0. The scalar helds with masses M+ and M_ are the messenger scalars. 



Figure 5: A Feynman diagram with the gaugino of mass m, the messenger fermion of mass 
M, and the messenger scalars of mass M±. 


the contribution from a U(l) gauge multiplet exchange. This U(l) gauge group is the toy- 
model version of the standard model gauge groups. The hat direction 0, 0 also has ±1 
charge under U(l), with D-hatness condition 0 = 0 . The result can be easily generalized 
to arbitrary gauge groups and messenger multiplets. We refer to the U(l) gauge coupling 
constant as g. The U(l) gauge multiplet acquires a mass m = 2g\(l)\. The task is to calculate 
the vacuum energy as a function of M, B, g and m. 

The Feynman diagrams are shown in Figs. In all our calculations we expand the 
amplitudes in terms of MB/M"^ and keep only the leading non-trivial terms of 0{MB)‘^. 
We need MB/M'^ < 1 to avoid a color- or charge-breaking vacuum, and this expansion is 
known to be a good approximation for the supersymmetry breaking mass squared for the 
hat direction unless MB is very close to m- 

We start with diagrams (Al) and (A2) in Fig. 


(Al) = ig^ry 


f yp 

f yk 1 

1 

1 (27r)4j 

(27r)"^p2 — 



(A.l) 


26 









Figure 6: Feynman diagrams with vacuum polarization due to (Cl) messenger scalar loops, 
(C2) “seagull” diagram with messenger scalars, and (C3) messenger fermions. 


(A2) 


f d^p 

r d^k 

1 

1 {271)^ J 

{2nY 

_{k^ - Ml){p‘^ - Ml)_ 


(A.2) 


Since we know that the sum of all diagrams vanishes in the supersymmetric limit MB —> 
0, we subtract the corresponding amplitude in the supersymmetric limit from each diagrams. 
The diagrams (Al) and (A2) give after the subtraction: 


(A1),+ (A2), 



+ x(l 


d^P r d^k 

{2ny J {2nY p"^ — w? 

1 

x)p^ - M2 + (1 - 2x)MBf 


{MB 0) 


(A.3) 


Here and hereafter, the subscript s refers to the subtraction of amplitudes in the supersym¬ 
metric limit. 

We expand the integrand in powers of MB/M^. The linear terms in MB vanish upon x 
integration, and we are left with the following expression, to 0((MB)^) in the integrand. 


(A1),+(A2)^ = ig^ f dx 
^ 0 


d^p 


d^k {4-l){l-2xf{MBf 

(27r)"^p2 — {k^ -|- x{l — x)p^ — M'^Y 


(A.4) 


We follow the same strategy as above to compute the contribution from the diagram (B) 
in Fig. |] containing the messenger fermions. 


(B). 




d'^k 


2{k ■ p) 


{p + kf-Ml 


{2tiY J (27r)"^ — M^) 

+ {Ml Ml) - 2{Ml MY 


/ dx 


d'^p f d'^k xp^ 


{2ti)^ J (27r)"^ — m? 

1 


(F + x{l - x)p‘^ - M2 - a;MB)2 


+ {MB -MB) - 2{MB 



27 






















= 


1 f d^p f d^k 


dx 


p 




( 27 r)^ J ( 27 r)"^ {k‘^ + x{l — x)p‘^ — 


+ 0{B^). 


(A.5) 


Finally the diagrams (Cl), (C2), (C3) with the gauge boson loop in Fig. ^ We hrst 
calculate the vacuum polarization diagrams of messenger helds. Note that the contribution 
of messenger fermions (C3) is the same as the one in the supersymmetric limit, and hence 
cancels after the subtraction. The scalar loop gives 


r d'^k 

(Cl). = g^J-^,{2k + pr{2k + pr 

^ ^ +{Ml 


_k^ - Mlik + pf - Ml 
and the “seagull” diagram gives 




2{Ml 




(A.6) 


(C2), =-29^/ 


d^k 

w? 


k^-Mi 


+ {Ml Ml) - 2{Ml 


M^ 


(A.7) 


Their sum is 


{Cl)^ + {C2)^ = g\-g'^y+p^pn 


d^k B 


dz 


Q{1 - 2z)‘^{MBf 


{2'k)'^ Jo {kJ + z{l — z)p^ — M'^y 


+ 0{B^). (A.8) 


Now including the gauge boson loop, the total contribution of the vacuum polarization 
diagrams is 


r d^k 

f d^p f 

1 (27r)4j 

{2n)^Jo 


dz- 


p 


{i-2zy 


(C), = t9g\MBy^ 


Adding all diagrams, we obtain —iVeff{m‘^) = (A)^ + (B)^ + (C)^ up to 0{MB)‘^, 


(A.9) 


Vejf{m^) = -l2g‘^{MBf f dx 

0 



r d^p 

f d^k 

1 

1 

to 

to 

1 

to 

H 

CO 

_1 

{2ti)^J 

to 

1 

to 

{P -I- x(l — x)p‘^ — M^)^ 


(A.IO) 

After a Wick rotation in the k and p spaces, the d'^k integration can be carried out trivially. 
Note also that the denominator is symmetric under the interchange of a: <->■ (1 — a:). We 
therefore symmetrize the x integration by substituting the polynomial in the numerator 
P{x) = (1 — 2a;)^ — 2a;^ by \{P{x) + P(1 — x)) = —x(l — x), hnding 


Veff{m^) = 


2 


-g\MB) 

12SM Jo 


dx / p^d{p^) 


x{l — x)p^ 


(p2 + m?){x{l — x)p^ + M'^y 


(A.ll) 


The p^ integral is logarithmically divergent. Fortunately, the divergent piece is nP indepen¬ 
dent, and hence is the renormalization of the cosmological constant. We subtract Veif(O) 


























from the above expression and redefine it as Veff. The hnal integral is convergent for any 

G [0, oo): 


Veff{m^) = 


1287r4 


pi poo 

dx p^d{p^ 
h Jo 


a;(l — x) 


(p2 + m'^){x{l — x)p^ + M2)2 


(A.12) 


The p integration can be computed using the following tricks. First, change the integration 
variable to = a;(l — x)p^. Then the integration can be done in the standard way, and 
we obtain 




g\MBf 

B 

l/z^ — x(l — x) + x(l — x) ln[x(l — x)z‘^] 

1287r4 

1 UjJU 

Jo 

(l/z^ — x(l — x)y 


(A.13) 


Here and below, we use = m?/M"^. 

We can further perform the x integration using dilogarithms. Using the roots of the 
denominator a = (^1 — — A //2 and 1 — a. 


Veff {a) 


B [(x - a)(a; - 1 + a) + a;(l - x) In (gfzfy) 

1287r4 Jo (x-a)2(x-l + a)2 


(A.14) 


After the hnal integral is carried out we are left with an expression for the effective potential 
as a function of a: 


Veff (a) 


g\MBf I 

I 

1 -2a(l - 
(1 - 2a) 


ln(a(l — a)) 
(l-2a)2 


-a) 


In^(a) 



a) — Li 2 (ci) + Li 2 (l 



(A.15) 


The form of the effective potential is shown in Fig. as a function of z as a solid line. 

The expression is manifestly real for z^ > 4. In the limit z —*• cxo, a ~ Ijz^ and the 
potential behaves as 


Veff = 


rsj 


g‘^{MBf 
647r^ 

g\MBf / 
1287r4 \ 


ihiV)-ln(^") + ^ + 0 
2 0 


In 


4g2 


eM2 



(A.16) 


This asymptotic form of the effective potential is also shown in Fig. as a dotted line. 

In the case < 4, a is complex, a = 1/2 + — 1. We can make the effective 

potential manifestly real, using the following dilog relations: 

Li 2 (a) + ^ln^(a) = -LU (-) - ivr ln(a) + ^ - C(2), (A.17) 

2 \aj 2 

Li 2 (l - a) + ^ ln^(l - a) = -LU + ztt ln(l - a) + ^ - C(2). (A.18) 

^ \ X Oj / Zi 
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Figure 7: A plot of the effective potential Eq. ( |A.15|) for 2 ; = 2g\(j)\/M < 30, in the unit of 
g‘^{MBY/ (1287r^). The solid line is the exact result, and the dotted line shows the asymptotic 
form Eq. (|A.16| ) valid for large 2 ;. 


We find 


V.if(z^ < 4 ) 


g\MBf 




^TT + arctan 




(A.19) 


Note that Li 2 ( 2 ;) — Li 2 (^*) is pure imaginary. In the above form, it is simple to take the limit 
> 0, and we obtain 


Veff = 


9\MBf 


■Hz^) + 


ln{z^) + "-^-z-"-^-z)+ 0{z^ \n{z^)) 




(A.20) 


= 4 — 


a \2 /MB 


4:71/ \ M 


+ |.^|^) + 0(4>% 


(A.21) 


The approximate form Eq. (|A.2CI|) truncated at 0{z‘^) is shown in Fig. ^ as a dotted line 
together with the exact form Eq. ( A.15|) (or equivalently, Eq. ( A.19|) ) as a solid line. From 
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Figure 8: A plot of the effective potential Eq. ( |A.15| ) (or equivalently, Eq. (|A.19| )) for the 
small held amplitude, z = 2g\(j)\/M < 1, in the unit of g‘^{MBY The solid line is 
the exact result, and the dotted line shows the approximate form Eq. (|A.20|) valid for small 


the last expression ( |A.21|) with a = g'^/ dvr, we can read off the mass of the hat direction. For 
messengers in 5 + 5* representation, we multiply the hnal result by a group theory factor 
r“T^tr(T“T^) = ^Cf where the trace is taken over the messenger helds and Cf is the second 


order Casimir for the hat direction. We obtain = 2Cf 
that in Ref. IHI 


- 

47r / 




which agrees with 


B Effective Potential and Wave-function Renormaliza¬ 
tion 

When one computes an ehective potential, one can determine the location of the minimum. It 
is well-known that one also needs to evaluate the wave-function renormalization 
in order to discuss the time evolution of the scalar held in general. Fortunately, such a 
calculation is not necessary in our case. 

Let us recall the simple fact that the ehective potential in our case is at 2-loop order: 
V ~ (a/dvr)^. Since it is a hat direction in the supersymmetric limit, this is the only term 
in the potential. The equation of motion in the hat space is 

Here we dropped the friction term 3H(j), but the essence of the following discussions does 
not depend on such simplifying assumptions. 
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Because V is of order (a/dTr)^, the motion is suppressed by a power in the coupling 
constant. Note that Z\(l)) is at most order (a/dTr). By factoring out the coupling constant 
factors, 



(B.2) 


C(0), 

(B.3) 

we hnd 

It is convenient to rescale the time variable t by 


(B.4) 

a 


(B.5) 

and we hnd 

<9^ / a \ C'(0) ( dcpy 

(9r^ ) Z{(j)) ytlr j 


(B.6) 


It is clear that the leading terms in the equation of motion are given by d‘^(j)/dT‘^ +u'(0) = 0, 
and all dependences on the wave function renormalization occur only at higher orders in 
perturbation theory. Therefore, the calculation of the effective potential is enough for our 
purpose, and we do not need 2’(0). 


C Time Evolution of the Flat Direction 


The evolution of the flat direction 0 is interesting in the LEGM models. Once the amplitude 
is dominated by the gauge-mediated piece, the potential is approximately proportional to 
(In and the dilution of the coherent oscillation occurs much slower than in the parabolic 

potential case. In this Appendix, we investigate the evolution of the flat direction by using 
the virial theorem. 

The virial theorem tells us that. 


2{K) 




(C.l) 


where K = (f>*(f> is the kinetic energy. In our case, V ~ Vo(ln |0p/(5')^)^ with (S) ~ 
3 X 10"^ GeV, and 

For In \(1)\‘^/{S)^ 3> 1, the energy density of the held is potential dominated. 
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The field equation is 


0 + - 0. (C.3) 

Multiplying it by 0* and using the energy density £ = K + V, 

Since the energy density is dominated by the potential term, we can write down the approx¬ 
imate evolution equation of the amplitude |0|, 


d\(l>\ orrui 

and hence 

\(j){t)\R{t)^ ~ constant 

This formula is valid when In |0p/(5)^ 3> 1. 


(C.5) 

(C.6) 


D Estimation of the Primordial Baryon Asymmetry 

In this Appendix, we justify Linde’s formula, Eq. ( b-6|) in our case. In fact, the validity 
of the formula depends on the nature of the operator O and the time evolution of the fiat 
direction 0. We clarify the reason why the formula is valid in the cases of our interest. We 
start from the equation of motion for the baryon number density, Eq. (|5.5|), 


riB + SHub = i 


d(j)' 



(D.l) 


It is useful to rewrite the equation in terms of baryon-to-entropy ratio Yb = tib/ s, to find 


V 1 



(D.2) 


where we used the relation sR^ = constant. Assuming a vanishing initial value YB{to) = 0, 
we obtain 

1 fdO do \ 

A crucial question is whether the t integral is dominated a.t t ^ to or t ^ oo. 

In the following analysis, we assume the Universe to be radiation dominated when the 
field begins to roll down the potential, R oc Another assumption is that the baryon- 
number violating operator O can be treated as a small perturbation to the evolution of the 
field 0. 
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In the case 0o ~ 0eq, or in the general hidden sector case, we have O cx 0^. The evolntion 
of 0 is essentially determined by Tny24>'^ by assumption and hence 0 oc oc On 

the other hand, s oc R~^ oc The integrand in Eq. (|D.3[) therefore behaves as and 

hence it is dominated by f ~ fg- By putting them together. 


Tb(oo) = 



dO 


00 * 


f-3/2 

+- 3/2 


1 8ml,2 


lm(0^) 


2m2/2lm(0^) 


(D.4) 


which essentially justihes Eq. (|^). 

In the other case of interest, 0o^0eq, both the behavior of the operator and time- 
evolution are completely different as discussed in detail in Section ^ We have (9 oc 0^, while 
0 oc R~^ oc Then the integrand behaves as which is unexpectedly the same as 

in the previous case. By putting them together, we obtain 



1 8EoIm(0^) 2EoIm(0^) 

2 ” 


(D.6) 


which again essentially justihes Eq. ( ^.61) . 

As it is clear from above the derivations, Eq. (^.6|) is not necessarily valid if the integral 
is dominated at f ~ oo rather than f ~ fg- We have not seen an explicit discussion on this 
point in the literature. 


E Affleck-Dine Baryogenesis in Hidden Sector Scenario 

In this appendix, we discuss the Affleck-Dine baryogenesis based on the scenario with SUSY 
breaking in the hidden sector. In this case, the gravitino mass is much larger than in the 
LEGM case, and all the scalar helds also have the SUSY breaking masses of the order of 
the gravitino massj^ In particular, the potential for the hat direction is always given by the 
supergravity contribution, which is essentially parabolic with a curvature of the order of the 
gravitino mass. 


(E.1) 

with m 3/2 ~ 1 TeV. Due to this fact, the evolution of the hat direction is much simpler than 
in the LEGM case. 

*In this Appendix, we denote all the soft SUSY breaking masses for the scalar fields by TO3/2 for simplicity. 
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Even if the gravitino mass is about 1 TeV, Eq. (|5.6| ) is still valid since the baryon number 
is generated when 0 starts to oscillate. With the baryon number violating operator (|5.4|) , 
we obtain 


ne Ho 


'”^3/2 


^^ 3/2 100 

mT 


sin 400- 


(E.2) 


and hence 


ns 
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-l/4_ 


100 


Ho 


'^■3/2 


1/2 ,^ 7/2 


sin 46*0- 


(E.3) 


If there is no entropy production after this stage, the above formula gives us the resulting 
baryon-to-entropy ratio. 

If there is entropy production, the primordial baryon number is diluted. The primary 
source of the entropy is the decay of the flat direction. Here, the potential for the flat 
direction is always parabolic, and 0 starts to oscillate when T = Tq ~ 
as discussed in Section |^. Then, by using the relation |0pT“^ = const., the background 
temperature at the 0 decay is given by 


^dec ~ g* 



(E.4) 


where 0dec is the amplitude of the flat direction when it decays. Furthermore, the reheating 
temperature due to the decay of 0, Tr, is given by 


~ S'* ^'^^\/"^3/20dec 

Then, the dilution factor is given by 


Ti 

TL ~ 


(E.5) 


(E.6) 


Usually, 0 
decay rate 


decays when the expansion rate of the Universe, if, becomes comparable to the 
of 0, T,^. In Ref. |]TI| , T,^ is estimated as 




rsj 



and hence ii ~ T^ results in 





(E.7) 


(E.8) 
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Combining the above results, we obtain 


Ub 

s 


a.y/» / Af. y-'* 

TT/ V”^3/2/ 
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sin 46*0 = 
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/ »^3/2 \ 

U TeV; 


- 1/6 

sin 46*0- 


(E.9) 

Another potential source of entropy is the Polonyi held related to the SUSY breaking, 
or the moduli helds in the string theory, which also have masses of order 7713/2 • The critical 
difference between the hat direction and the Polonyi held z is the formula of their decay 
width; since the Polonyi held couples to particles in the observable sector only through 
interactions suppressed by the gravitational scale, its decay width P^ is much smaller than 
the width of (j). As discussed in Section ^ P^ is estimated as 


Nch ^1/2 

Att M2 ■ 
* 


(E.IO) 


Even with this decay rate, we can apply an argument similar to the case of the entropy 
production due to 0 decay; Eqs. ( |E.2| ) and (|E.3| ) are still valid, and we also obtain equations 
similar to Eqs. (|E.4|) - ( |E.6|) where 0’s are replace by z’s. The remainder is to evaluate the 
amplitude of z at its decay time, Zdec, by using the relevant formula for P^. By solving the 
equation M ~ P^ with H ~ m 3 / 2 ^/M*, we obtain 


^dec ~ 


^1/2 
M, ■ 


(E.ll) 


Then, assuming the initial amplitude of z to be zq ~ M*, the dilution factor is given by 
D ~ M^/m^/ 2 , and hence 


ub im3/2lm(0^) 

- ~ J-y -::- 

S M 2 


rsj 


/I0oiy 
V M J {Mt ) 


sin 46*0- 


(E.12) 


Thus, the baryon-to-entropy ratio may be larger than ~ 10“^° even after the decay of the 
Polonyi held. 

However the reheating temperature after the decay of 2 ; is likely to be too low. By using 
Eq. ((E.IOI) , the reheating temperature is estimated as 

Tr ~ ~ IMeV x ■ (E.13) 

Thus, if the gravitino mass is heavier than about 10 TeV, the Polonyi held may decay 
before the big-bang nucleosynthesis (BBN), and the scenario which gives Eq. (|E.12|) may be 
viable. [] However, for a favorable range of the gravitino mass ( 7773/2 ^ ITeV), the reheating 

Hn fact, even if m3/2 ^ lOTeV, there may still be a problem since the lightest superparticle produced by 
the decay of the Polonyi held may overdose the Universe |^, To solve this difficulty, we may have to 
accept a much larger gravitino mass, or a scenario in which the lightest superparticle in the MSSM sector is 
unstable. 
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temperature is less than 100 keV which is lower than the temperature where the big-bang 
nucleosynthesis (BBN) starts. This means that the decay of z signihcantly affects the results 
of the standard BBN scenario. In this case, we need some mechanism to reduce the energy 
density of the Polonyi held. A thermal inhation [|^ is an interesting candidate for it. The 
baryon-to-entropy ratio in this case is discussed in Section ^ By using the fact that the 
ratio of ub to pz is constant in time, we obtain 


Ub 

s 


rsj 


Pz 


X m 


-1 

3/2 



4 

sin 46*0. 


(E.14) 


Thus, once we hx the ratio pz/s after the late inhation, the baryon-to-entropy ratio can be 
estimated. 
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